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. When neutrinos are produced from long lived particles like pions, kaons 

qq ■ and muons, the life-time of the source particle affects the flavour conversion 

probability formula. Experiments like LSND which use muon decay neutrinos 
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are two orders of magnitude more sensitive to lower values of mass square dif- 
ference compared to other experiments where the sources are pions or kaons. 

The standard formula for neutrino oscillations is derived with the assumption [|TJ that the 
uncertainity in the initial position of the neutrino is small compared to the distance between 
$_i 1 the production and the detection sites of the neutrinos . We shall show here that in a 
covariant treatment, the effective spread of the netutrino wavefunction in space is given by 
[a x + cr) where a x is the initial spatial uncertainity of the production point of the neutrinos 
and r is the lifetime of the sources (eg. pions, kaons , muons etc) whose decay produce the 
neutrinos. We find that the effect of long lifetime sources is to exponentially suppress the 
oscillation term in the conversion probability fromula. The flavor conversion probability as 
a function of distance , for relativistic neutrinos produced from long lived resonances, turns 
out to be 

P(u e -> X) = - sm 2 26 (1 - cos{ }exp{-{ ) 2 }) (1) 
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where Am 2 is the mass square difference in eV 2 ,L where Am 2 is the mass square difference 
in eV 2 ,X is the detector distance in m (km), r is the lifetime of the source particle in the lab 
frame in m (km) and E is the energy in MeV (GeV). When cr is comparable to the detector 
distance spatial oscillations of the conversions probability is not seen. However the neutrino 
mass squared difference that can be probed is lower for longer lived sources. We find that for 
the LSND experiment where the neutrinos are produced from muons (cr M = 658.65m), when 
fitted with the covariant oscillation formula (|l|) gives a bound on Am 2 which is two orders 
of magnitude lower than the corresponding bounds obtained from other similar experiments 
but where the neutrino sources are pions or kaons (cr n = 7.8m). 

We start with an initial wave function of a general Gaussian form 

*?(x,t) = 1 exp{ i(P a (x - Xi ) - E a (t - U)) - - (2) 

The initial spread of the wavefunction in space around the mean initial position Xj is denoted 
by a x and the spread in the time direction around the mean initial time of production tj is 
denoted by a t . The magnitudes of a x and a t depends on how the state is prepared. The 
time evolved wavefunction can be obtained from the initial wavefunction (0) by the linear 
superposition principle , which can be written formally as 

ty a (x f ,t f ) = J dxdt K(x f -x,t f -t) ^(x,t) (3) 

where K (xf —x,tf — t) is the probability amplitude of a particle initially located at (x, t) to 
be detected at another spacetime point (x/, tf). The transition amplitude K(xf — x,tf — t) 
for a free particle is given by the expression ( |2| [TJ) 

K(x f -x,t f -t;m a ) = (—) (-^) i^i(im a s) (4) 

where s = ((tf — t) 2 — (xj — x) 2 ) 1 ^ 2 is the invariant spacetime interval propagated by the vi 
mass eigenstate. If this interval is large (s >> m~ Y ) and time-like ( (tf — t) > (xf — x) ) 
then we can use the asymptotic expansion of the Bessel function 

/ 2 \ 1/2 

KAims) ~ ( ) exp{— ims} (5) 

\irimsj 
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to obtain from the expression for the propagation amplitude at large time-like separation 

K(x f -x,t f -t;m a ) = ( — ^ ffl ) 1/2 exp{-im a ^J (t f - tf - (x f - x) 2 } 

2ni^J (tf — t) — (xf — x) 

(6) 

This expression for the free particle propagator is valid for bosons. For fermions their is a 
extra factor (i jd + m a ) operating on the l.h.s of the expressions (f|) and (|6]). This factor only 
changes the normalisation and the expressions for the conversion probability is identical for 
bosons and fermions. The time evolved wave-function can be evaluated by substituting the 
expression @ for the propagator K(xf — x,tf — t) in the expression ([|) for 
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^a{ x f ~ x i)tf ~ U] m a) = N J dx dt 



a/4 



2 - 



xex P mx,t)- - ^-^} (7) 



where the phase factor in the exponential is of the form 



t) = -m ay J{ff - tf - (x f - x) 2 + P a (x - x % ) - E a (t - U) (8) 

and the constant coefficients have been clubbed together as the factor N which normalises 
the wave-function. We perform the integrations in (0) by the method of stationary phases 
||. The integral is approximated by the integrand along the trajectory where the phase is 
an extremum. The extremum of the phase (|j) is given by (d<&/dt) = => (m a (tf — t)/((tf — 
t) 2 - ( Xf - x) 2 ) 1 ' 2 ) = E a and (d$/dx) = (m a (x f - x)/{{t f - t) 2 - (x f - x) 2 ) 1 ' 2 ) = P a . 
One can solve for t and x using these two equations and substitute in the integrand of (^) 
to get the expression for the integral in the stationary phase approximation , which turns 
out to be 

((Xf - Xi) - (t f - tj)#-) 2 
* a (x/ - x^ tf - U) = N exp {-iE a (t f - U) + iP a (xf - Xi) +a 2tPa\2\ a J 

V x t \ E a ) ) 

(9) 

where again we have clubbed together all the constants as the normalisation coefficient N'. 
Denoting X — Xf — x i: T — tf — and v a = (P a / E a ) which can be identifired by distance of 



propagation, time of propagation and the group velocity of the particle respectively we can 
write the expression for the time evloved wave-function compactly as 

* a (X, T) = N'exp {—iE a T + iP a X - ( * ~ ^ ' } (10) 

4(0-2 + ^2 a 2) 

The state vector of a mass eigenstate can be in the mass basis \v a > in the form 

\m a ;X,T>=t> a (X,T)\u a > (11) 

On the other hand the state of a weak interaction state say v e is expressed as a linear 
combination 

\v e ]X,T >= ]>> a (X,T)K >< z/> e > (12) 

a 

The probability amplitute of a neutrino to be produced at (x i? ij) as a v e and to be detetected 
at (xf,tf) as another weak eigenstate is given by 

A{p e ^v^X,T) =< z/> e ;X,T>=^ a (X,T) < p„\u a >< v a \v e > (13) 

a 

and the corresponding probability is given by 

P{X, T; v e -> ^) = I < z/> e ; X,T > | 2 

= ^^ a (X,T)^*(X,T) < ^|i/ a >< i/ a |i/ e > <v e \v h xv h \v il > (14) 
0,6 

Restricting ourselves to mixing between two flavour generations, and denoting the mixing 
matrix elments as < u a \u e >= — < u b \u^ >= sinO and < v^v^ >=< u h \u e >= cos9, the 
expression (|i"4|) for the probability reduces to the form 

P{v e -> u„X,T) = sin 2 9 cos 2 9 {V a (X, T)%(X, T) + %(X, T)%(X, T)} 

- sm 2 cos 2 # {^ a (X, T)^*(X, T) + ^(X, T)* 6 (X, T)} (15) 

This expression expresses the probability of a z/ e produced at (x i; £,) to be converted to a z/ M 
at some other spactime point In actual experiments only the distance X between 
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the source and the detector is known but the arrival time of the neutrinos is not measured. 
The expression for the conversion probability as a function of only X is obtained by taking 
the time average of ([H]). The normalisation is chosen so that *ff a (X,T) represents a single 
particle wave-function and therefore, 

/+oo 
dT* a (X,T)** a (X,T) = l (16) 
-oo 

The non-trivial contribution to the time averaging comes from the interference term of ([15]) 

/+oo r+oo 
dTV a (X, T)** b (X, T) = Re / dT exp {-t(E a - E b )T + i(P a - P b )X} 
-oo J — oo 

This integral can be evaluated by the completing the squares, to give 

Re [ + °° dT^ a (X,T)%(X,T) = cos{((E a - E b )^f^- - {P a - P b ))X} 

7-oo {vl + v£) 



x ■ ,.2n ((^ - E >) ° - T^r )} (18) 

where 



K + v b ) 

We express the energy, momentum and the masses in terms of their averages E = (E a + 
E b )/2,P = (P a + P b )/2,m = (m„ + m b )/2 and differences AE = (E a - E b ),AP = (P a - 
P b ), Am = (m a — m b ). In terms of these variables (|18|) reduces to the form 

V2 //rnArn _ A 



, - cos ( (EAE - PAP)-) e~ A , (20) 

{v* + vi)y 2 p 

/ • _S H O" , . , AP^ . r> X 2 . . 

A = ( ) y(^) + ("p-) a 4i5 ( 21 ) 

where we have retained the terms to the first order in AE/E, AP/P and Am/m. The energy 
and momenta of the neutrinos are determined by the energy-mommentum conservation 
laws at the production vertex. Since niether the energy nor the momenta of the remaining 



outgoing particles are measured, one cannot fix either the energy or the momenta of the 
different neutrino mass eigenstates in the linear combination state. So one cannot assume 
either E a = Ef, or P a = P^. Only the mass shell relation E 2 = P 2 + m 2 a for each mass 
eigenstate v a ( P||). In terms of the average and difference , the mass shell conditions 



imply, 



A 777 ^ 

EAE - PAP = (22) 



Using the relation ([22|) in fl2~i~D we see that the interference term reduces to the form 

\/2 . Am 2 i 4 .„„. 

cos ( — —X) e~ A (23) 



(v 2 a + v 2 y/ 2 v 2P 

we see that the oscillation length is L osc = (4irP/Am 2 ) for both relativistic as well as non- 
relativistic particles. The significant difference in the formula comes from the suppression 
factor A fl2~T|). For intereference term to be non-zero, both the terms of the suppression 
factor must be small. The second term in ( pip is << 1 as long as X « L co h with the 
coherence length for relativistic neutrinos given by , 

E 2 

L coh = AV2(a x + v<r t ) (^) (24) 

where v = (v 2 + v 2 )/(v a + v^) ~ 1. This expression for the coherence length differs from 
the expression derived in the standard wave-packet treatments U by the presence of the a t 
term. In all accelerator neutrino oscillation experiments the criterion X « L co h is satisfied 
, so the effect of this term is negligible. The first term of (|2T|) plays a more significant role. 
For relativistic neutrinos, the suppresion factor (|21|) reduces to the form, 



A = (^ ^L±^))3 (25 ) 

In the accelerator experiments o x - the spread in the beam of primary particles is of the 
order of a few cm. The neutrinos are produced from the secondary decays of pions and kaons 
produced in the primary collision. The uncertainity in time of production of such neutrinos 
at is then given by the lifetime r of the pions, kaons or muons which are the source of 
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neutrinos for that particular experiment. In such situations a x << a t and [o x + va t ) — vr 
and the expression fl2"5|) reduces to 

A - i^f. (26) 

Using these results we see that the expression for the time average of the conversion 

probabilty, for relativistic neutrinos produced from long lived resonances is given by 

^/ . 1 9 „ . .Am 2 X. ,Am 2 T,, , , . 

P(v„ -> z/ e ; X) = - S m 2 2fl (1 - co S (^^) exp - {^j^f ) (27) 

where Am 2 is the mass square difference ,X is the detector distance , r is the lifetime of the 
source particle in the lab frame and E is the neutrino energy. In practice one averages over 
the energy flux n(E) of the neutrinos and the average probability which is fitted with the 
experimental number to obtain the allowed regions of Am 2 and sin 2 28 is given formally by 
the expression, 

The limits on the values of Am 2 and sin 2 29 obtained by fitting the results of different 
experiments with the covariant oscillation formula (P7|) are listed in Table I and plotted 
in Fig 1, Fig 2 and Fig 3. The regions of the parameter space allowed by the covariant 
wavepacket formula and by the standard formula of the LSND (muon source) experiment 
are shown in Fig.l. and the Karmen (pion source) experiment are shown in Fig 2. Since 
the pion decay length is of the same order as the experimental baseline, the improvement in 
sensitivity to Am 2 is marginal compared to LSND. The combined result of all experiments 
listed in Table I plotted using the covariant oscillation formula is shown in Fig 3. 

The dependence of the oscillation term on the source lifetime can be understood as 
follows. When the time uncertainity is large the uncertainity in energy becomes small and 
the wave-functions are then eigenstates of energy. In the time smearing the overlap of two 
different energy eigenstates dissapears and therefore the oscillation term vanishes. 

Acknowledgments I thank Terry Goldman, Walter Grimus, Harry Lipkin and Eligio Lisi 
for their comments on the earlier papers P|§. 



TABLES 

TABLE I. The asymptotic limits on Am 2 and sin 2 29 from different experments according to 
the oscillation formula (27). r is the lifetime of the neutrino source in the lab frame, < E u > is 
the average v energysource in the lab frame, < E v > is the average v energy, X is the detector 
distance and P is the experimental value of the conversion probability. 



Experiment (Source) 


t (m) 


< E v > (MeV) 


X (m) 


P 


Am 2 (eV 2 ) 


sin 2 26 


LSND ( M ) | 


658.6 


30 


30 


(0.16 - 0.47) x 10" 2 


(1.0 - 1.5) x 10" 3 


0.003 - 0.009 


LSND (tt) jn| 


17 


130 


30 


(0.26 ±0.15) x 10^2 


0.4-0.8 


0.002 - 0.0082 


Karmen (tt) [| 


7.8 


29.8 


17.5 


< 0.3 x 10" 2 


< 0.08 


< 0.6 x IO" 2 


E776 (tt) |[ 


578 


5 x 10 3 


10 3 


< 0.15 x 10~ 2 


< 0.1 


< 0.3 x 10~ 2 


CCFR( K) U 


5.41 x 10 3 


140 x 10 3 


1.4 x 10 3 


< 0.9 x 10" 3 


< 0.6 


< 0.18 x 10" 3 


Bugcy(U,Pu) 


3 x 10 10 


5 


95 


< 0.75 x 10- 1 


< io- 9 


< 0.15 
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FIG. 1. Lsnd decay at rest experiment j|] allowed regions with the standard formula 
(between the dotted curves)and the covariant oscillation formula (continuous curves). 
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FIG. 2. Karmen tt decay at rest experiment Q allowed regions with the standard formula 
(dotted line) and the covariant oscillation formulas (dashed curve) . 
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FIG. 3. The region between the continuous lines is allowed by the LSND fi experiment p| .The 



region between the dotted lines is allowed by the LSND tt experiment [11] . Region ruled out by 
E776 j7j is above the dashed-dotted curve and by Karmen Q is above dashed curve. The region 
above the top-most dashed curved is ruled out by CCFR ||. 
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